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ABSTRACT 

We present a study of the core collapse supernova (CCSN) rate in nuclei A and Bl, 
of the luminous infrared galaxy (LIRG) Arp 299, based on ~ 11 years of Very Large 
Array (VLA) monitoring of their radio emission at 8.4 GHz. Significant variations in 
the nuclear radio flux density can be used to identify the CCSN activity in the absence 
of high-resolution very long baseline interferometry (VLBI) observations. In the case 
of the Bl-nucleus, the small variations in its measured diffuse (synchrotron plus free- 
free) radio emission are below the fluxes expected from radio supernovae (RSNe), thus 
making it well-suited to detect RSNe through flux density variability. In fact, we find 
strong evidence for at least three RSNe this way, which results in a lower limit for the 
CCSN rate (^ SN ) of > 0.28±0.16yr _1 . This value agrees within the uncertainties with 
the infrared (IR) luminosity based SN rate estimate, and with previously reported 
radio estimates. In the A-nucleus, we did not detect any significant variability and 
found a SN detection threshold luminosity of « 3.1 x 10 28 ergs _1 Hz -1 , allowing only 
the detection of the most luminous RSNe known. Our method is basically blind to 
normal CCSN explosions occurring within the A-nucleus, which result in too small 
variations in the nuclear flux density, remaining diluted by the strong diffuse emis- 
sion of the nucleus itself. Additionally, we have attempted to find near-infrared (NIR) 
counterparts for the earlier reported RSNe in the Arp 299 nucleus A, by comparing 
NIR adaptive optics images from the Gemini-N telescope with contemporaneous ob- 
servations from the European VLBI Network (EVN). However, we were not able to 
detect NIR counterparts for the reported radio SNe within the innermost regions of 
nucleus A. While our NIR observations were sensitive to typical CCSNe at ^300 mas 
(or 70 pc projected distance) from the centre of the nucleus A, suffering from extinc- 
tion up to Ay ~ 15 mag, they were not sensitive to such highly obscured SNe within 
the innermost nuclear regions where most of the EVN sources were detected. 

Key words: galaxies: individual: NGC 3690:IC 694 - galaxies: starburst - stars: 
supernovae: general 



1 INTRODUCTION 

Galaxies with infrared luminosities > 10 11 Lq, i.e., lumi- 
nous infrared galaxies (LIRGs) show enhanced activity as- 
sociated ^^h_&_^axbw^_wc]J_a£a.ri active galactic nucleus 
(AGN) (|Sanders &: Mirabell [l996). These galaxies contain 
a high concentration of dense molecular gas in their cir- 
cumnuclear regions; consequently a very active star for- 
mation process is ensured, resulting in a large number of 
core collapse supernovae (CCSNe). However, the detection 



of the SNe occurring in the nuclear regions of LIRGs is 
in most cases not feasible at optical wavelengths due to 
severe dust extinction (see SjSJ. Observations in the near- 
infrared (NIR) K-hanA where the extinction is strongly re- 
duced (Ak ~ 0.1 x Ay) can be used fo r their success- 
ful d etection (e.g.. iMattila fc Meikiel l200ll : iMaiolino et all 
2002). However, the tight concentration of the star forma- 
tion within the i nnermost nuclear regions of LIRGs (e.g., 
ISoifer et all l200ll ) means that the detection of such SNe 
even at NIR wavelengths in natural seeing conditions is of- 
ten not feasible. Recent SN search programmes benefiting 
from the high spatial resolution offered by adaptive optics 
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(AO) observations at NIR wavelengths have already discov- 
ered several obscured SNe close to LIRG nuclei. For exam- 
le, SN 2004ip with a likely high extinction Av > 5 mag 



Mattila et al.ll2007l ) was detected at a projected distance of 
1.4arcsec, or 500 pc, from the rf-band nucleus of the nearby 
LIRG, IRAS 18293—3413. Its core-collap se nature was con- 
firme d by a VLA detection at 8.4 GHz |Perez- Torres et al.l 
120071). Two circumn uclear SNe, SN2004iq and SN2008cs 
( Kankare et al. ! 2008) were detected in the circumnuclear re- 
gions of the nearby LIRG IRAS 17138-1017. SN 2008cs was 
shown to suffer from a record high extinction of Ay ~ 18. 
Its core-collapse nature was agai n confirmed by a radio de - 
tection at 22.4 GHz using VLA l| Perez- Torres et al] liooli) . 
SN 2010cu detected at 0.4arcsec (1 80 pc) from the jf- band 
nucleus of the nearby LIRG IC883 dRvder et alJl201Ch . has 
also a likely high extinction. More recently, another proba- 
ble SN (PSN J13203538+3408222) was detected at 0.8 arcsec 
(360 pc) from the centre of t he if-band nucleus of the same 
galaxy (|Kankare et aljfeoill ). 

Whilst assuming a realistic initial mass function (IMF), 
the detection of new CCSNe can yield information on the 
star formation rate (SFR ) of the host galaxy for the progen- 
itor stars (Mnit > 8M0. ISmartt et al.| [2009). provided that 
the number of detections over time allows a reasonable sta- 
tistical estimate. Many of the CCSNe expected to occur in 
LIRGs will be radio emitters, and their evolution is deter- 
mined by the circumstellar medium (CSM) in which they 
are embedded. Despite the expectation of a large number of 
CCSN events taking place in LIRGs, only a few tens of them 
have been detected at radio frequencies (free of extinction), 
either as young supernova remnants (SNRs) or as RSNe in 
the nuclear regions of e.g . , Arp 220 and Arp 299 (|Parra et al.l 
l2007l : IPerez-Torres et afll2009l . and references therein). The 
difficulty in their detection is due to the presence of a strong 
host galaxy background emission. Very long baseline inter- 
ferometry (VLBI) observations provide high angular resolu- 
tion and are able to resolve out the diffuse background emis- 
sion, therefore allowing the detection of RSNe and SNRs in 
the innermost nuclear regions of nearby LIRGs. 

The VLA does not provide us with the needed reso- 
lution to resolve SNe in the innermost nuclear regions of 
LIRGs. However, it is possible to use VLA observations to 
search for variations in the total nuclear radio flux density, 
which can be used to probe the RSN activity in LIRGs. 
Based on the variability of the flux density, we propose a 
method to estimate CCSN rates (^sn) in nearby LIRGs by 
exploiting the wealth of archival VLA data which are avail- 
able for a number of LIRGs. Archival VLA data for Arp 299 
at 8.4 GHz, covers a period of time suitable to test the pro- 
posed method and to allow a meaningful comparison with 
other f S N predictions. 
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Figure 1. Deep image of the 8.46 GHz radio continuum of 
Arp 299 resulting from combining (in the «?>-plane) the archival 
VLA observations (11 epochs) which have an image rms < 
150 /^Jy (see Table gj. Cores A, Bl, C, and C are indi- 
cated, as well as the probable background source D (see Ta- 
ble [TJ. The locations of previous optical and NIR SNe (see 
Table [3j are indicated by five-pointed stars. The contours are 
(3,3 \/3,9,...)x45/iJy beam" 1 , the off- source rms flux density per 
beam. The restoring beam size is 0.54x 0.52 arcsec with a position 
angle of 80° . 



Table 1. Positions of the Arp 299 cores measured from the image 
shown in Figure [T] and their uncertainties. The position of the 
probable background source D is also given. 



Core 


Peak Position 


A[a,6] 
(mas) 


a(J2000) 


<5(J2000) 


A 


11 28 33.626 


58 33 46.65 


0.2 


Bl 


11 28 30.987 


58 33 40.80 


1.6 


c 


11 28 31.332 


58 33 50.00 


4.5 


c 


11 28 30.648 


58 33 49.30 


5.8 


B 


11 28 33.013 


58 33 36.55 


5.6 



(NGC 3690) and C' (in the overlapping region). Hereafter we 
will refer to the NGC 3690-nucleus as Bl, since this is the 
brightest core in that galaxy. In the same field of view a prob- 
able background source named D is also found, which to our 
knowledge has not been detected in any radio spectral line 
observation made towar ds the system (see e.g.. ICasoli et al] 
Il999l : lAalto et all 119971 ) . indicating that its redshift differs 
from that of Arp 299. In fact, source D wa s first detected at 



radio wavelengths b vlHuang et al.l dl990t) and identified as 
a probable SN, but lUlvestadl l)2009l ) has found that D has 
characteristics of a Palomar-Green (PG) quasar. The posi- 
tions of the cores in Arp 299 and the probable background 
source D, as measured from the deep combined VLA image 
at 8.46 GHz (see Figured]), are listed in Table[j] 



At a luminosity distance of 44.8 Mpc ( Fixsen et al.l 



2 PREVIOUS RADIO AND NIR STUDIES OF 
ARP 299 

Arp 299 is a system composed of the pair of galaxies 
IC694 (east component) and NGC 3690 (west component), 
whose interaction is induc ing a powerful star formation 
l|Alonso-Herrero et alJ I2000L and references therein) . Four 
separate cores can be distinguished in VLA radio images of 
Arp 299 (see Figure [TJ: A (IC694), B (= Bl + B2) plus C 



Table 2. Radio based SN rate estimates for nuclei A and Bl. 



i^sn (yr 1 ) References Notes 

A Bl 



0.65 0.12 Alonso-Hc rrero et al. ('2000') 

0.5-1.0 0.1-0.2 Neff et al. (2004) 

1 The estimates are based on VLA measurements. 

2 No errors have been reported for the given rates. 
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1 19961 ) for H = 73 kms^ 1 Mpc~\ Arp 299 has an infrared 
(IR) luminosity (L IR = L[8-1000/xm]) « 6.7 x 10 11 L 
(as resulting from scalin g the value for the same wave- 
length interval given by ISanders et all [2003, to the dis- 
tance we adopted here), which makes it one of the most 
lumin ous LIRGs of the local Universe. ICharmandaris et al.l 
|2002i ) found that IC 694 emits approximately 40% of the 
total IR luminosity of Arp 299, whereas NGC3690 ac- 
counts for approximately 20%. Making use of these esti- 
mates together with the empirical relation between CCSN 
rate and I/[8-1000Miri] for starburst galaxies obtained by 
iMattila fc Meikld (200ll ): 



yr- 1 



2.7 x 10" 



we have a total CCSN rate of 1.8 yr~ for the whole sys- 
tem, and we expect SN rates of ~ 0.7yr _1 and ~ 0.4 yr -1 
for A and Bl, respectively. These IR luminosity based SN 
rates are similar for the nucleus A, but higher for nucleus 
Bl, than those determined in previous studies (see Table 
[2J, which bear some dependency on the radio spectral be- 
haviour of the cores. We notice though, that these IR lumi- 
nosity based SN rate estimates assume no AGN contribution 
to the IR luminosity. However, we now know that there is 
likel y also some AGN contribu tion to the IR luminosity of 
Bl (|Alonso-Herrero et al.l [20091 1. and therefore the true SN 
rate for this nucleus could be lower than what has been es- 
timated above. 

Over the past 20 years, Arp 299 has been the target 
of many SN search campaigns which can also provide 
important insight into its dominating energy source 
(AGN and/or starburst). In that p eriod a number of 
SNe hav e been detected a t optic al llTreffers et al.l 1 19931 : 
iLi et all Il998l: iQiu et all 1 19991: iNewton et all l201Ch 
or NIR Jvan Buren et al.1 1 19941 : IMattila et"aU 120051 : 
IMattila fc Kankard l2010h wavelengths (see Table [3]) in 
the circumnuclcar regions. All of these SNe have been 
spectroscopically classified as core-collapse events of either 
type II or lb, except SN 1992bu, whose type is unknown. 
The detection of seven SNe over a period of less than 20 
years, already indicates a very high SN rate in the circum- 
nuclear regions of Arp 299. The relatively high fraction of 
stripped envelope SNe (types Ib/IIb) and their distribution 
with respect to the ho st galaxy properties are discussed in 
lAnderson et alJ (|201ll ). 

Ver y high r e soluti on radio studies b y iNeff et al.l 



20041): luivestadl J200! 



iPerez-Torres et al.l |2009] ) and 
iPerez- Torres et all (|2010T ) have shown that at high angu- 
lar resolution the brightest cores of the Arp 299 merger (A 
and Bl), consist of a wealth of compact components within 
the innermos t regions of these n uclei. 

In 2008, IPerez-Torres et ail started an observing cam- 
paign with the EVN to monitor the nuclear region of 
Arp 299-A. The extended radio emission of Arp 299-A is 
resolved out with the VLBI beam, thus revealing the ex- 
istence of a rich cluster of 26 compact radio emitting 
sources within the 150 pc diameter innermost region. The 
high brightness temperatures observed for the compact 
sources are indicative of a non-thermal origin for the ob- 
served radio emission, imp lying that most of the sou rces are 
youn g RSNe and SNRs (IPerez-Torres et all l2009l l. More- 
over, IPerez-Torres et aU l|2010l ) have recently detected the 



Table 3. Optical and NIR Supernova history of Arp 299. The 
references covering discovery, type determination and position 
are given. 



SN Type 
designation 



Position 
5 (J2000) 



1992bu (NIR) 1 
1993G (optical) 2 
1998T (optical) 5 
1999D (optical) 7 
2005U (NIR) 10 
2010O (optical) 12 
2010P (NIR) 14 



IIL 3 

lb 6 

II 8 

lib 11 

lb 13 

Ib/IIb 15 



11 28 31.46, 58 33 38.3 1 
11 28 33.43, 58 33 31.0 4 



11 28 33.16, 
11 28 28.42, 
11 28 33.22, 
11 28 33.86, 
11 28 31.38. 



58 33 43.7 6 
58 33 38.7 9 
58 33 42. 5 10 
58 33 51.6 12 
58 33 49.3 14 



van Buren et al.lll994l 
Treffers et al.lll993l 



Tsvetkovlll994l 
Forti et al Jjl993l 



Li et alj|l998f 



Yamaoka et al 


Il998l 


Oiu et al.l 


1999 




Jha et al.l 


1999 



Prospcri 1999 



Mattila et aT1|2005 
Modiaz et"aIll200E 




iLeonard & Cenkoll2005l 


Newton et al|l201( 




Mattila et al.ll201C 


Mattila & Kankardl20ld 


Rvder et al.ll20ld 



AGN in the A-nucleus, whose characteristics correspond to 
a low luminosity AGN (LLAGN). In the western nucleus of 
Arp 299, Bl, on ly five VLBI components have been detected 
jUlvestadll2009l) . 

However, Arp 299 has been observe d with VLBI tech- 
niques only since 2003 |Neff et al.ll2004l ). and thus the ap- 
pearance of new sources in its nuclear regions before 2003 
has gone unnoticed. We have therefore carried out a study 
using the available data in the VLA archive aiming at mea- 
suring the ^sn in the Arp 299 nuclei through the study of 
their flux density variability. Our method is particularly ef- 
ficient in the Bl-nucleus, where the radio surface brightness 
is lower, hence allowing us to detect small variations in the 
overall nuclear flux density, which can be related to SN ac- 
tivity. In fact, this kind of study has been successfully ap- 
plied to M82, where a strong radio flare on top of the diffuse 
emission of the host galaxy, was identified with SN 2008iz 
jMarchili et al.ll2010l ). 



3 VLA ARCHIVAL DATA 
REDUCTION AND ANALYSIS 

Arp 299 has been observed repeatedly with the VLA over the 
last decades in X-band (8.46 GHz). However, archival data 
lack homogeneity in frequency and resolution and do not 
cover an adequate time period at regular intervals, making 
the estimates of Vsn more challenging. 

The X-band VLA observations in A-configuration, pro- 
vide us with a resolution of a few tenths of an arcsecond, 
enough to easily distinguish the emission from the different 
cores in Arp 299, but not to resolve each nucleus. In fact, 
as we have mentioned in jJU t he nuclei A and Bl both con- 
sist of several compact sources (I Perez- Torres et alj|2009l . and 
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Table 4. Properties of the VLA images and estimated flux densities at X-band (y ~ 8.46 GHz), Stot, for the radio light curves 
of nuclei A and Bl (see Figure f2(a)[ |. The peak intensities, Speak, f° r each nucleus are also given. The flux density and peak 
intensity uncertainties have been computed by adding in quadrature the rms noise in the map (sixth column), plus a conservative 
5% uncertainty in the point source calibration. We have used a common uv-r&nge of 23.6 to 261.9 kA for all the epochs, and 
0.78 X 0.73arcsec at 19° as the convolving beam. 



Date VLA Programme Weighting rms ^peak S iot ^pcak S foi 











config. 






(/^Jy/beam) 


(mjy/bcam) 




(mjy) 




(mjy/beam) 


(mjy) 


1 


1990- 


-02- 


-24 


A 


AH396 


natural 


148.4 


64 


,74±3 


.24 


83 


,66±4. 


19 


7.52±0.40 


9.89±0.52 


2 


1990- 


-03- 


•04 


A 


AH396 


natural 


89.3 


62 


,77±3 


.14 


79 


,46±3. 


97 


7.55±0.39 


9.60±0.49 


3 


1991- 


-07- 


•05 


A 


AS333 


natural 


93.3 


62 


,26±3 


.11 


81 


,87±4. 


09 


8.09±0.41 


11.23±0.57 


1 


1993- 


-01- 


-28 


BnA 


AS333 


natural 


73.4 


62 


,32±3 


.12 


81 


.71±4. 


09 


7.91±0.40 


12.29±0.62 


5 


1993- 


-05- 


•07 


B 


AS333 


uniform 


294.8 


61 


,56±3 


.09 


78 


,94±3. 


96 


7.95±0.50 


11.34±0.64 


6 


1994- 


-05- 


•16 


BnA 


AY064 


natural 


38.6 


64 


,20±3 


.21 


82 


,12±4. 


11 


7.81±0.39 


10.39±0.52 


7 


1999- 


-10- 


-28 


B 


AS568 


uniform 


70.8 


56 


,25±2 


.81 


77 


21±3. 


86 


7.57±0.38 


9.92±0.50 


8 


2001- 


-03- 


•19 


B 


AN095 


uniform 


37.4 


63 


,35±3 


.17 


83 


,03±4. 


15 


8.38±0.42 


12.48±0.62 


9 


2002- 


-04- 


29 


A 


AN103 


natural 


207.7 


63 


08±3 


.16 


82 


,78±4. 


14 


7.14±0.41 


10.46±0.56 


10 


2002- 


-09- 


•06 


B 


TYP04 


uniform 


332.1 


60 


,79±3. 


.06 


76 


,37±3. 


83 


7.35±0.50 


10.25±0.61 


11 


2003- 


-10- 


•17 


BnA 


AS779 


uniform 


36.9 


66 


,94±3 


.35 


82 


,54±4. 


13 


7.97±0.40 


11.48±0.57 


12 


2004- 


-11- 


-02 


A 


AC749 


natural 


31.5 


65 


,85±3 


.29 


84 


.37±4. 


22 


7.80±0.39 


10.90±0.55 


13 


2005- 


-02- 


-08 


BnA 


AW641 


uniform 


103.5 


65 


,46±3 


.27 


84. 


,58±4. 


23 


8.19±0.42 


12.23±0.62 


14 


2006- 


-04- 


-15 


A 


AC749 


natural 


131.6 


58 


,83±2 


.94 


77 


39±3. 


87 


15.47±0.78 


17.28±0.87 



references therein). Therefore, care must be taken while in- 
terpreting flux density measurements at VLA resolutions, 
which do not resolve out each core into their individual 
sources. 

Our aim is to further investigate the available archival 
data to extract information on the SN activity in Arp 299, 
despite the difficulties both in extracting and interpreting 
the information from VLA images. We have analysed high 
resolution, X-band VLA data from 1990 to 2006 (see Table 
[4| to study the flux density variations of the Arp 299 nu- 
clei. The VLA data were reduced following standard proce- 
dures with the NRAO Astronomical Image Processing Sys- 
tem (AIPS). The intensity scale was set by observations of 
3C48 in the case of project TYP04, and 3C286 for all the 
other projects. At X-band, these sources have an adopted 
flux d ensity around 3.2 and 5.2 Jy, respectively (|Baars et al.l 
1 19771 ). 

M ost of the s t udied epo chs have been p reviously anal- 
ysed bv lNeff et al] (|2004n and lUlvestadl (|2009T l. The main dif- 
ference with the previous studies, is that we have estimated 
the flux densities in a consistent way by extracting them 
from maps created with matched baselines (in wavelengths) 
and with the same convolving beam (0.78x0.73 arcsec at 
19°). We chose a common arrange of 23.6 to 261.9 kA for 
all our maps in order to properly compare among different 
observing epochs, thus avoiding systematic effects. Such re- 
striction in the wy-range implies a suppression of structures 
larger than «8. 7 arcsec and smaller than «0. 8 arcsec. More- 
over, to match the resolutions of the different observations, 
we applied natural weighting to those images with better 
resolution (those made in the VLA-A configuration) and 
uniform weighting to images with worse resolution (those 
made in the VLA-B configuration). For the hybrid BnA ob- 
servations, we applied a uniform weighting to images having 
an intrinsically worse resolution than the chosen convolving 
beam, and a natural weighting to images having a better 
resolution than this. 

We ran a series of phase-only self-calibration iterations 
in each epoch to achieve a better signal to noise ratio. All 



the images were centred at ll' l 28 m 33 s , +58°33'43", and the 
AIPS task PBCOR was applied in each one of the final im- 
ages to correct for the primary beam attenuation. 

We estimated the flux densities of the nuclei in each im- 
age by enclosing the emission above 7a in a box and fitting a 
single Gaussian, using the AIPS task IMFIT. In each epoch, 
such box had an area of « 500 x 500 pc 2 for the A-nucleus, 
and « 400 x 400 pc 2 for the Bl-nucleus, thus covering the 
circumnuclear regions in each case, where we expect the SNe 
and the possible AGN activity to take place. The angular 
size of these regions is such that fitting a single Gaussian is 
suitable. 

Table[4]contains the observing parameters of each epoch 
and the estimated values obtained always from the last im- 
age resulting from the self-calibration process, but corrected 
by the beam attenuation. 



4 TRACING THE NUCLEAR ACTIVITY IN 
ARP 299 THROUGH FLUX DENSITY 
VARIATIONS 

Figure [2 ( a) | shows the light curves of nuclei A and Bl, scaled 
to the same flux density level (the one of Bl) to allow com- 
parison. We note that the Bl-nucleus exhibits conspicuous 
flux density variations. In contrast, the A-nucleus does not 
show systematic variations and remains relatively constant 
throughout the period of time comprised in our study, tak- 
ing into account the observational uncertainties. We have 
also attested that the use of different phase calibrators has 
no influence on the flux density variations we see in Bl. 
Furthermore, as we show in Figure |2(a)| the flux density 
variations of nucleus Bl are not correlated with those of nu- 
cleus A, and we have checked that the same is also true with 
the other cores present in the region, i.e., C, C' and D. 

To quantify the flux density variations we see in nucleus 
Bl, we have made an iterative selection of the lowest data 
points to obtain a baseline value for Bl, and compared each 
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epoch's flux density with respect to the obtained baseline 
value. 

As a starting point, we define an initial baseline value 
for the Bl-nucleus flux densities given by (S) ± 8S, where 
(S) is the mean of the 14 Stot values and SS is the standard 
error of the mean. We then determined the deviation of the 
flux density in each epoch (i) from the baseline value, with 



<Stot, 



(8) 



8S 



The same was done for nucleus A. In the next step we ex- 
cluded the epochs having dSi > 3 and proceeded to calculate 
a new baseline value. We followed this iterative process until 
no more epochs could be removed by another iteration, thus 
leading to a robust baseline value. 

After the third iteration, we find that the flux den- 
sities of the Bl-nucleus on epochs 1, 2, 6, 7, 9 and 10, 
are representative of a quiescent state. Thus, considering 
the fluxes of these epochs, we obtain a final baseline value 
of (5 B1 ) ± SS B1 = 10.08 ± 0.22 mJy for Bl. For nucleus 
A, only epochs 13 and 12 were found to be above three 
times the standard error of the mean after the second iter- 
ation, and we have consequently obtained a baseline value 
of (S A ) ± 8S A = 80.59 ± 1.16 mJy for A. Taking the ratio 
of the final A and Bl nuclei baselines values, we obtain a 
scaling factor for A of ~ 0.125, that allows comparison with 



Bl (see Figure 2(a) I 



To attest the significance in the variability of each epoch 
in each nucleus, we calculated again the deviation of each 
epoch but this time using the standard deviation, of the 
final baseline for each nucleus: a A = 2.57 mJy and <r B i = 
0.34 mJy. We find that the flux densities of all the epochs 
in the A-nucleus show no significant deviations from the 
mean value (< 3<r A ). For the Bl nucleus, epochs 3, 4, 5, 
8, 11, 13 and 14 (following the numeration from Table [4]) 
are all above 3cr B i, and of these, epochs 4, 8, 13 and 14 are 
also above 5a B i- Epoch 12, not being part of the baseline, 
has a variation slightly below 3cr B i, and hence it lacks of 
significance according to our detection threshold. Epochs 3 
and 5 appear to be connected with the variation present in 
epoch 4, so we consider all these three epochs as having the 
same origin. Epoch 11 has less significant variability on its 
own and seems to share a similar behaviour with the epoch 



11 in A (see Figures 2(b) and 2(c) I, so its significance is 
further reduced due to a possible systematic effect in that 
epoch. 



4.1 Flux density variability of the A-nucleus 

Tracing the SN activity of nucleus A through its radio flux 
density variations, is not feasible. This is because the flux 
density of the core A at these angular scales is not domi- 
nated by the individual SNe, but by the diffuse emission of 
the strong starburst therein, thus resulting in the relative 



quiescence of the total flux density (see Figure 2(b) I 



We note that typical RSN at the distance of Arp 299 
have flux densities comparable to the given uncertainties in 
the nuclear flux density (see column 7 of Table [4]), and thus 
to the obtained standard deviation, <r A , in our variability 
test. Hence, it was not possible to detect SNe in the A- 
nucleus using this method. Furthermore, the minimum flux 
density for an event to be robustly detected with 5<r A above 



the mean, would be ~ 12.85 mJy. At the distance of Arp 299, 
such flux density implies L ~ 3.1 x 10 28 erg s" 1 Hz" 1 at 
the peak, which is uncommon for normal CCSNe (see Fig- 
ure [3]). With our method, we are able to detect only very 
bright SN events (such as powerful type Iln or 98bw-like 
events), which do not seem to be common in the A-nucleus. 
Indeed, the lumin o sities for the SN candidates found by 
I Perez- Torres et al.l l|2009l ). range between w 3 x 10 26 and 
« 2 x 10 27 ergs" 1 Hz" 1 (typical of type lib, IIP and IIL 
SNe), and are thus well below our detection limit. Thus, in 
the absence of VLBI observations, the only feasible alter- 
native to detect new SNe in nucleus A is by means of NIR 
obs ervations (see §5§. 

iPerez- Torres et all (2010) identified a compact source 
as the core of a LLAGN in Arp299-A. Its flux density at 
5 GHz has remained below ~ 1 mJy in o bservations carried 
out i n 8 April 2008, 5 Dece mber 2008 dPerez-Torres et ail 
I2009T ) and in 7-8 June 2009 (| Perez- Torres et alj|20ld ). In 
fact, typical LLAGNs display at most flux density varia- 
tions ~ 0.8 times the average flux density (see Figure 2 from 
iNagar et al]|2002h . Even if the LLAGN doubles its radio lu- 
minosity at a given time, we would not be able to detect such 
a burst in nucleus A since our flux density uncertainties are 
higher than 2mJy in all the epochs (see Table [4]). 



4.2 The nature of the radio variability of nucleus 
Bl 

As mentioned above, epochs 4, 8, 13 and 14 show conspic- 



uous variability above 5<r B i (see Figure 2(c) I. Following the 



same reasoning as for the A-nucleus, we have also explored 
the possibility of explaining the flux density variations in the 
Bl-nucleus being due to AGN outbursts. The presence of a 
LLAGN in the Bl-nucleus is well supported by th e detection 
of H2O megamaser emission (jTarchi et al .1 [20T1I. and refer- 
ences therein) and hard X-ray comp onent (e.g.. iBallo et al.l 
|2004| ; iGonzalez-Martfn et all 120091 ) . Further evidence for 
a LLAGN in Bl comes from the s tudy of mid-infrared 
(MIR ) high-excitation emission lines (|Alonso-Herrero et al.l 
2009). He nce, o ne of the five VLBI components reported by 
Ulvestad (2009) (without considering that one correspond- 
ing to the 2005 event) is a candidate for being the LLAGN. 
Each one of these sources has a flux density 0.3 mJy at 
both 8.4 and 2.3 GHz, typical for LLAGN. In the event that 
one of these sources is the AGN, and even if it were to dis- 
play the highest inter- year variability th at has ever been 
observed from LLAGN IjNagar et al.ll2002i ). the correspond- 
ing burst will have a flux density <C 5a B i, thus being below 
our detection limit. Therefore, we find no evidence for any 
of the detected events to be linked to AGN activity. 

We consider next the possibility of explaining the nu- 
clear variability in 1993, 2001 and 2005 (epochs 3-5, 8 and 
13-14) as resulting from supernova explosions. Note that the 
scarcity of data available makes it challenging to properly 
determine the physical parameters that describe the vari- 
ability we have detected. Nevertheless, we can identify cer- 
tain properties in each event by examining their rise and fall 
in the light curve (Figures 2(a) and 2(c) I. Moreover, if we 



plot the peak spectral luminosity versus the time after the 
ex plosion date to reach the peak (as was done for SN2000ft 
bv lAlberdi et alj|2006h . we will find that the events consid- 
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Figure 2. |(a)| Arp299-A light curve (solid line with circular markers) and -Bl (solid line with squared markers) at X-band. The labelling 
numbers correspond to those of Table[4] The Arp299-A flux densities have been scaled to the Bl-nucleus baseline value (dotted line) for 
comparison. The scaling factor used is ~ 0.125. The epoch zero days corresponds to 1 January 1990. Figures | (b) | and | (c) | show the actual 
flux measurements for nuclei A and Bl, respectively. In both Figures, the solid line represents the baseline value (80.59 zt 1.16 mjy for 
nucleus A, and 10.08 ± 0.22 mjy for nucleus Bl) and the dotted lines denote the +3o" and —3a levels, a being the standard deviation: 
2.57 mjy for A and 0.34 mjy for Bl. We note that all the epochs in A are well contained within the ±3<ta region and essentially all the 
epochs represent the baseline level, whereas in Bl a number of epochs show significant variation > 3cr B i- 
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Table 5. Properties of the variability events Bl-nucleus found 
during a ~llyear period of VLA observations at 8.46 GHz. 



Event 


'Speak 


-^pcak 


^pcak 




(mjy) 


(10 27 ergs -1 Hz" 1 ) 


(years) 


1993 


> 2.21±0.66 


> 5.31±1.58 


<2 


2001 


> 2.40±0.66 


> 5.76±1.58 


<1.3 


2005 


> 7.20±0.90 


> 17.29±2.16 


>1.5 
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Figure 3. Peak spectral radio luminosity of different types of CC- 
SNe versus the time to reach the peak after the explosion date 
(updated from lAlberdi et al.l I2006T ) . We have added a factor of 
^/(5GHz) multiplying tpcak that allows us to plot I/ pea k a t any 
frequency. Labels for some historic RSNe are added for compar- 
ison with our 1993, 2001, and 2005 events (asterisks). Note that 
the asterisks refer to Z/ pea k 

lower limit values and t peak upper 
limit values for all events, except for 2005 where t pca ]j is a lower 
limit. The dotted line represents the minimum luminosity a SN 
should have to be detected at 5oa in nucleus A through flux den- 
sity variability. The dashed line represents such detection limit 
(5<t B i) for nucleus Bl. 



ered here, have typical values for moderately luminous type 
Ib/c and II SNe (see Figure |3j. 

A couple of assumptions must be made. First of all, 
we consider the peak flux density of the events (S pca k) as 
the difference between the flux density at the peak epoch 



(as seen in Figures 2(a) and 2(c) I, and the characteristic 
baseline value of the Bl-nucleus (10.08 ± 0.22 mjy), which 
we consider as the flux density of the nucleus in its quiescent 
state. Secondly, we take half the time between the rise and 
fall in the light curve as the time it takes for each of the 
putative SNe to reach their peak after the explosion date 
(tpeak)- Since each event is poorly sampled we can, as a 
coarse approach, consider their peak flux densities only as 
lower limits, and their time to reach the peak as upper limits. 
We summarize these two assumptions for each putative SN 
event in Table [5] 

From Table O we see that the peak flux densities of the 
Bl-nucleus corresponding to the 1993 and 2001 events, are 
above 2.2 mjy. At the distance of Arp 299 these flux densi- 
ties correspond to luminosities of the order of a few times 



10 27 erg s -1 Hz -1 , which are typical for type Ib/c or mod- 
erately luminous type II SNe, as we can also see in Figure 
The 2005 event is more luminous than the others and 
similar to the most luminous type II RSNe observed, e.g., 
SN 1986J and SN 1988Z. In the following, we discuss each 
event separately: 

• 1993 event (epochs 3-5). This event exhibits an ex- 
tremely slow rise and relatively fast fall. Typical RSNe do 
not display such behaviour. If a single SN is responsible 
for the increase in flux, it must have reached its peak be- 
tween epochs 3 and 4 (thus having a flux density peak 
higher than that of epoch 4). A change of the circum- 
stellar medium (CSM) profile could explain the changing 
slope in its optically thin phase. In this case, we assume 
that the putative SN has exploded some time after epoch 
2 and has reached its peak some time between epochs 3 
and 4, therefore being characterised by £ pca k < 2 years and 
Lpcak > 5.3 x 10 27 ergs -1 Hz -1 . These values for L pca k and 
tpcak are expected from moderately luminous type Ib/c or 
type II SNe, as we see in Figure [3] 

• 2001 event (epoch 8). This event is characterised by 
i P c a k < 1.3 years (obtained by taking half of the time be- 
tween epochs 7 and 9) and L poa k > 5.7 x 10 27 ergs -1 Hz -1 . 
This corresponds most likely to a type Ib/c or a moderately 
luminous type II SN (following the same criteria as before) . 

• 2005 event (epochs 13th-14th). This event corresponds 
to the appearance of a new compact source (SN) in the inner- 
most nu clear region of B l, as detected with VLBI observa- 
tions bv lUlvestadl(|2009D . The VLA observations we present 
here, cover just the beginning of the rise of this event, there- 
fore allowing us to solely estimate lower limits to its peak 
luminosity and time to reach the peak (see Table [5}. Yet, 
we can still place this event among the most luminous RSNe 
known to date (Figure [5}. 

Unfortunately, Arp 299 was not observed with the VLA 
at X-band between 1994 and 1999 and we cannot search 
for new SN explosions in the Bl-nucleus during that pe- 
riod. Hence, we divide our analysis in two periods: the first 
one from February 1990 to May 1994 (hereafter PI), and 
the second from October 1999 to April 2006 (P2). This is 
equivalent to one period of 4.2 yr and another of 6.5 yr. 

Event 1993 in PI, as well as events 2001 and 2005 in 
P2, can all be explained by RSNe. We note that the RSN 
type of each event is poorly determined due to the coarse 
and uneven sampling of the light curve. Only a more fre- 
quent sampling could provide the necessary information to 
characterise each eve nt, as was done for SN2008iz in M82 
iMarchili et al.l (|2010l ). We could solely state that due to the 
substantial difference in lifetime of type Ib/c and bright type 
II SNe, our events are likely to be type II SNe rather than 
type Ib/c. Nevertheless, the mere identification of events to 
have a CCSN origin, allows us to estimate the CCSN rate. 



4.3 The CCSN rate in nucleus Bl 

We now want to obtain an estimate for the ^ S n over the 
~ 11 years (P1+P2) comprised in our study. In order to do 
this, we assume that the N SN events detected in Bl are 
independent and occur at a constant rate over time, thus 
obeying Poisson statistics. Consequently, the uncertainty on 
this number is a =-\/~N. 
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In the case of PI, with one RSN (1993 event), we have 
^sn=i ~ 0.24 ± 0.24 yr -1 . For P2, where there are two de- 
tections (2001 and 2005 events), we have ^sn = 2 ~ 0.31 ± 
0.22 yr -1 . We then calculate a total ^sn for the ~11 years 
comprised in our study by taking the average of the PI and 
P2 estimates, which yields a final z^sn=3 ~ 0.28 ±0.16yr _1 . 



Our variability test leads us to a ^ S n estimate which, 
within the uncertainties, agrees with the IR luminosity based 
SN rate estimate, as well as with previous radio-based esti- 
mates (see Table[5]). We remark, however, that our estimate 
can only represent a lower limit due to our detection thresh- 
old. Event 1993 is only slightly above such limit, and we see 
from Figure[3]that fainter (as well as more rapidly evolving) 
SNe than the one responsible for the 1993 event, are deemed 
to remain as non-detections through this variability test. We 
probably miss most t ype Ib/c and all type II P SNe, which 
repre sent about 36% (Li et al.ll2oTiT ) and 59% (jSmartt et all 
120091 ). respectively, of all CCSNe in a volume- limited sample 
(see Figure [3| . Thus, the total i'sn should be significantly 
higher than our estimate, unless the IMF is a top-heavy 
one favouring the production of the most massive stars, 
and the surrounding environment is dense enough, so that a 
large fraction of the SN explosions will turn into brig ht ra- 
dio sources. For instance. IChevalier fc Franssonl (|200lh have 
shown that the high molecular cloud densities in the nu- 
clear regions of galaxies undergoing bursts of star formation 
(e.g., M82), can account for the high radio luminosity of the 
SNRs therein. This is also the case for the luminous RSNe in 
Arp220, which iParra et all (20071 ') have found to be consis- 
tent with powerful type Iln SNe interacting with their dense 
CSM. In fact, the rate of appearance of the luminous RSNe 
in Arp 220 agre e s with the rate of 4±2 SNe per year found by 
lLonsdale et ail (20061 ), which implies a SFR high enough to 
power the entire Lir, of the galaxy. Therefore, if the gas den- 
sities surrounding the SNe in Arp299-Bl are high enough, 
our i^sn estimate could be close to the true CCSN rate. 
This result is not surprising if we conside r that already from 
the E VN observations of the A-nucleus, I Perez- Torres et al.l 
(2009) found that its u SK might be much higher than previ- 
ous estimates, thus suggesting that a top-heavy IMF could 
describe better the observations of the nuclear starburst in 
Arp 299-A. We might well be facing the same situation here 
for the Bl-nucleus. In fact, a recent study of the properties 
of the sample of optical/NIR detected SNe in the circumnu- 
clear regions of Arp 299, also indicates an I MF possibly bi- 
assed toward high mass stars in the system (Anderson et al.l 
H). 



However, the population of SNe within the innermost 
~ 200 pc of the nuclei of Arp 299 can be quite different from 
the one within ^1-4 k pc distances from the nuclei of Arp 299 
that lAnderson et al.l have studied, due to the difference in 
the amb ient conditions, which could lead to different IMFs 
(see e.g. IChevalier &: Franssonl l200ll ). Therefore, in the fu- 
ture the combined use of high spatial resolution radio and 
near-IR observations in the detection and study of the pop- 
ulation of SNe within the LIRG innermost nuclear regions 
is crucial in order to better constrain their star formation 
properties. 



4.4 Caveats and limitations of our variability test 

We have presented a variability test that covers a period of 
~ 11 years (PI + P2), with 14 data points. This corresponds 
to a typical sampling rate of 1.3 observations per year, and 
a typical time interval between observations of « 11 months 
(considering PI and P2 separately); although we acknowl- 
edge that the sampling is very uneven. 

The typical time interval between observations in our 
study makes the detection of type Ib/c SNe difficult, as they 
decay fast at radio wavelengths, having short lifetimes that 
range between a few days, up to a few months (see e.g., 
IWeiler et ai1l2002r i. It is thus impossible to resolve individ- 
ual type Ib/c SNe occurring within an 11 months period. 
However, it is not impossible to detect them, if we are for- 
tunate enough to have an observation at around the time 
of their luminosity peak. This is of course not very likely 
and our detection threshold (see Figure [3]) makes our test 
blind to practically all type Ib/c SNe (as we have discussed 
in H4.3|) . for the given sample. 

For type II SNe the situation is different. Although 
a coarse and uneven sampling does not allow an identi- 
fication of the SN subtype (IIP, IIL, lib, Iln), both the 
sensitivity and the frequency of our sampling, favour the 
detection of bright and long-lasting events. For instance, 
SN 2000ft in a LIRG host galaxy (NGC 7469) at 70 Mpc dis- 
tance, with a peak luminosity of almost 1 x 10 28 erg s _1 Hz -1 , 
was cle arly visible in VLA obse rvations over a period of 
4 years IjPerez- Torres et al. 2009). Furthermore, an event 
like SN 2000ft, could be covered by four of our observations 
(with a typical separation of ~ 11 months). In fact, having 
at least three observations covering the light curve, would 
be required for its minimum sampling. This is fulfilled in 
the case of the Bl-nucleus for long-lasting events. However, 
in the case of nucleus A the identification of RSNe is more 
challenging than in the case of Bl due to our current sam- 
pling characteristics, in addition to the sensitivity threshold 
which allows only the detection of very bright events. Ac- 
cording to IR luminosity based estimates, the SN rate of 
nucleus A (~ 0.7yr _1 ) is about twice that of Bl. A sam- 
pling rate of 1.3 observations per year is clearly not optimal 
to probe the RSNe activity in nucleus A, for which we could 
often have overlapping events. 



5 SEARCHING FOR NIR COUNTERPARTS 
OF RADIO SUPERNOVAE IN ARP 299-A 

NIR observations towards nuclei A and Bl were carried out 
with the Gemini-North Telescope in the A"-band (2.2 fim) 
under the program 'Altair Study of Supernovae in Lumi- 
nous Infrared Galaxies" (P.I.: S. Ryder). Our NIR observa- 
tions of Arp299-Bl will be reported elsewhere and here we 
concentrate on the analysis of the NIR data of nucleus A. 

In the previous section we have estimated the radio SN 
rate for the Bl-nucleus and pointed out that SN detections 
in the A-nucleus — unless extremely radio luminous SNe, like 
SN 1998bw — are feasible only by means of VLBI observa- 
tions, at radio wavelengths, or by NIR observations with 
high resolution. We then consider the EV N observations car- 
ried o ut in April and December 2008 by iPerez- Torres et al.l 
(2009), revealing young RSN candidates within the nucleus 
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Table 6. Log of the Gemini-North NIR observations of A- 
nuclcus. 



UT date 


Program 


Exposure Time (s) 


2008 Jun 20.3 


GN-2008A-Q-38 


8x30 


2010 May 04.8 


GN-2010A-Q-40 


8x30 



Table 7. Sources of uncertainty in the A-nucleus Gemini-N image 



wcs. 


Transformation 


RMS (arcsec) 


WCS of SDSS 


0.090 


SDSS to HST 


0.056 


HST to Gemini-N 


0.023 


Combined uncertainty 


0.1085 



A. Aiming at finding the NIR counterparts for these objects, 
we analysed our contemporaneous (June 2008) NIR adaptive 
optics (AO) Gemini-N observations. 

We used the Near-InfraRed Imager (NIRI) with the 
ALTAIR Laser Guide Star AO system (0.022 arcsec/pixel, 
FWHM~0.1 arcsec) for the two observations (see Table [6]) 
in this study. The obtained images were reduced using the 
NIRI package in IRAF V2.14. Before median combining the 
target images, we subtracted the inherent horizontal noise 
pattern of the NIRI images. 

The World Coordinate System (WCS) for the NIRI im- 
ages with a small field of view (23 x 23 arcsec) was obtained 
in an iterative way. We first calibrated the WCS of an SDSS 
i-band image using 28 point-like sources from the 2MASS 
catalogue. Having a total of 51 2MASS sources in the SDSS 
field, we rejected 23 which were either elongated, too faint, 
too close to another object to be blended, and/or located on 
corrupted pixels regions of the images. 

The SDSS i-band image was then transformed to an 
HST ACS J-band image (obtained as already calibrated 
from the HST science archive) using nine common sources 
between the two images. As a final step, we transformed the 
HST image to the images of the A-nucleus obtained with the 
Gemini-N telescope. For doing this, we used 20 sources. In 
all the steps we used a simple geometric transformation for 
the coordinates, which included shifts in x and y as well as 
the same scaling factors and rotation in x and y. The propa- 
gation of uncertainties between the different alignment steps 
yielded a total RMS of approximately 0.1 arcsec (see Table 

The centroid position (J2000) of nucleus A as mea- 
sured in the WCS calibrated Gemini-N image using IRAF 



is ll h 28 m 33!62, 58°33'46'.'6 (see deep image in Figure [4(a) I. 
The distance between the NIR A-nucleus and its radio coun- 
terpart (see Table [T]) is 0.103 arcsec, so that both positions 
are coincident within la. 



5.1 SN detection in NIR images 

We have used a difference imaging techni que based on the 
Optim al Image Subtraction (OIS) code (|Alard fe Luptonl 
1998), and already demonstra ted to work in SN detection 
in NIR AO images of LIRGs l|Mattila et all 120071 ). Having 



transformed the image from the epoch of interest (20 June 
2008) to our reference image from 2010 May 4, we then se- 
lected 19 regions around suitable bright objects for deriving 
a convolution kernel that is later used to match the point- 
spread function (PSF) of our reference image to the PSF of 
the other image before performing the image subtraction. To 
obtain a successful subtraction, we had to include also the 
nucleus itself in these regions, given its outstanding bright- 
ness compared to the other objects in the image. 

The FWHM of the A-nucleus is about 0.4 arcsec in our 
Gemini-N images. To obtain a satisfactory subtraction re- 
sult we used ~ 0.66 x 0.66 arcsec for the size of the regions 
involved in the kernel solution. This is large enough to cover 
well the diffraction limited core of the PSF (whereas the 
much wider PSF wings are hidden in the noise for most of 
the sources present in our images). Regions substantially 
smaller than this did not result in a smooth subtraction for 
the regions around the bright nucleus. As a result of this, the 
detection of SNe falling within the innermost ~ 300 mas ra- 
dius region is not feasible using these techniques (see Figure 



4(b)! 



EVN radio images bv iPerez-Torres et all fedOSl) have 
approximately 10 times better resolution than that of our 
NIR AO observations. Thus, the resolution of these NIR ob- 
servations is not high enough, such that we could expect to 
be able to detect objects as deep in the LIRG nucleus as 
in the EVN observations. In addition, SN detections in the 
NIR domain might also be hindered by the high extinction 
towards the SNe. However, not all the CCSNe produce radio 
emission strong enough to be detectable at the distance of 
the nearest LIRGs. Hence, we have investigated the limit- 
ing magnitude for our NIR SN detection ability at different 
distances from the centre of the nucleus A. 

5.2 SN detection threshold in the A-nucleus 

We have simulated SNe of different magnitudes at differ- 
ent distances from the centre of the nucleus A in our 2008 
June 20 image using the IRAF MKOBJECTS task. Unfor- 
tunately, there are no bright field stars in our Gemini-N 
images (FOV~ 23 x 23 arcsec) that can be used to model 
the complicated PSF. Therefore, for the fake supernovae a 
moffat profile with a beta parameter of 2.5 and FWHM of 
~ 0.15 arcsec was used. We found this to give a reasonable 
match with the diffraction limited part of the PSF whereas 
for the faint sources the much wider PSF wings disappear 
within the noise. 

After the matching and subtraction of the reference im- 
age for each simulation, we checked if the simulated source 
could be detected significantly (more than 3cr) above the 
subtraction residuals within the nuclear regions. The level 
of the noise introduced by the subtraction residuals was de- 
termined as the standard deviation of the aperture sums 
(using a small radius of 0.13 arcsec) for 49 apertures placed 
in a grid covering the innermost 0.66 x 0.66 arcsec nuclear re- 
gion. We obtained a photometric calibration for our Gemini- 
N images, with too small a FOV to include any 2MASS 
stars, making use of the HST photometry for several com- 
pact circumnuclear super st ar clusters (A1-A6) reported by 
lAlonso-Herrero et al.l (|2000l ) and converted to the standard 
if-band. We estimated our SN detection threshold by sim- 
ulating (one at a time) and recovering eight artificial SNe 
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Figure 4. Arp299-A WCS calibrated Gemini-N i^-band images (FOV= 15x15 arcsec), [(a)] Our reference image (June 2008), which we 
use to compare with contemporaneous VLBI images, A b l ow-u p of the inner nuclear region is also shown. The 'x' symbols indicate the 
positions of the EVN detected SNe by iPerez-Torres et al.l j2009l ) . [(b) | A subtraction of our Gcmini-N image obtained in May 2010 from 
the June 2008 epoch (shown as negative). Note that the use of the innermost part of the A- nucleus is needed for the PSF matching to 
obtain a smooth subtraction. However, because of this, we lose information about any possible SNe within a radius of ~ 70 pc (~ 300 mas) 
from the centre of the A-nucleus. SNe occurring outside the ~ 300 mas radius, arc detectable by our method with a limiting magnitude of 
~ 18.5, and outside 1 arcsec radius, to a limiting magnitude of ~ 20. Note that SN 2010O, discovered in January 2010 in the circumnuclear 
regions, is still well detected in the subtracted image. 



of a given magnitude, distributed at different locations, at a 
300 mas distance from the centre of the nucleus A (or ~ 70 pc 
projected distance). By repeating the simulation for a range 
of magnitudes, we found a 3<r limiting magnitude of ~18.5 
(estimating our calibration to be accurate to ±0.1 mag) for 
our SN detection. Inside this radius our method was not 
sensitive for the detection of highly obscured SNe. Note 
that the limiting magnitude we have found is similar to the 
magnitude of SN 20 04ip discovered wit h the VLT/NACO in 
IRAS 18293-3413 l|Mattila et al.ll2007n . 

We tested our SN detection ability also in the circum- 
nuclear regions where no strong image subtraction residu- 
als are present. For this the level of noise was first deter- 
mined as the standard deviation of aperture fluxes for 45 
apertures (0.13 arcsec radius, 0.20-0.26 arcsec sky annulus) 
placed along a circumference at 1 arcsec distance from the 
centre of the nucleus A. Again, by simulating and recovering 
12 artificial SNe of a given magnitude distributed at differ- 
ent locations along the same circumference, we found a 3cr 
limiting magnitude of ~ 20. 

From t he template light curve s for ordinary CCSNe 
provided bv lMattila fc Meikleil200ll . we adopt M(K) pcak = 
— 18.6 and a decline rate of Am(K) = 0.022 mag/day af- 
ter the peak. At the distance of Arp299, this corresponds 
to m(K) = 14.7 at the maximum light. We would there- 
fore be able to detect such a SN suffering from a visual 
extinction of up to Av ~ 15 (similar to the extinction found 
in our study for SN 2008cs in the circumnuclear regions of 
a LIRG) at 300 mas from the centre of nucleus A, still at 
three months from the maximum light. In the case of a 
slowly declining (SNe 1979C and 1998S-like) supernova, for 



which M(K) pea ) IL = —20.0, the apparent magnitude at the 
distance of Arp299 becomes m(K) = 13.3 at the maximum 
light. Thus the visual extinction towards such a SN could 
be as high as about 30 magnitudes and still detectable at 
300 mas from the centre of the nucleus A, three months af- 
ter the maximum. At a larger distance (1 arcsec or ~ 200 pc 
projected distance) from the centre of the nucleus, our obser- 
vations were sensitive to an ordinary CCSN with Ay ~ 30, 
or slowly declining CCSN with Av over 40. 

At radio frequencies there is no dust extinction, and a 
supernova will appear bright when the interaction between 
the SN ejecta and the CSM is strong enough, i.e., when 
the CSM is dense enough. However, at NIR wavelengths, all 
types of supernova events should be visible as long as the 
visual extinction is not too high. Therefore, supernovae in 
the less extinguished circumnuclear regions could be more 
easily detected at NIR wavelengths than at radio frequen- 
cies, while the opposite is true when dealing with objects 
within the innermost nuclear regions. 



6 DISCUSSION AND SUMMARY 

We have used VLA archival data to estimate the radio SN 
rate in the nearby LIRG Arp299. The wealth of data avail- 
able for Arp 299 in the VLA archive, makes this galaxy a 
unique target to perform such studies. 

We have used the variability of the nuclear radio flux 
density to probe the RSN activity in the absence of high- 
angular resolution VLBI observations. The radio flux densi- 
ties of the nuclei A and Bl within Arp 299, result from the 
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contribution of several compact sources located in the in- 
nermost nuclear regions, along with diffuse emission. In the 
case of nucleus A, we see that its flux density remains con- 
stant within the observational uncertainties throughout the 
time inte r val w e have studied here. This is consistent with 
lUlvestadl l|2009h claim that for this nucleus, the integrated 
flux density of the compact sources within the nucleus, is 
only about 20% of the total flux density seen at lower reso- 
lution. The rest is diffuse emission of the galaxy itself. 

Nevertheless, in the case of the Bl-nucleus, small vari- 
ations in flux density due to new supernovae can be distin- 
guished from the total nuclear flux density at VLA resolu- 
tion, or the SNe can even dominate it, as seen in the case 
of the 2005 event. In fact, we have been able to indirectly 
detect three radio luminous SNe within the Bl-nucleus, thus 
estimating a radio SN rate for Bl of i/gN > 0.28 ± 0.16 yr -1 . 
This estimate is in good agreement with the IR luminos- 
ity based SN rate estimate, and also with values obtained 
through other methods at radio wavelengths. The remark- 
able difference between our estimate and those obtained in 
other studies is that here, we have estimated the SN rate by 
directly identifying the RSN activity within the Bl-nucleus, 
which is not possible at other wavelengths. Note that the 
fact that any putative AGN is extremely weak (as sug- 
gested by VLBI observations) , allow us to interpret flux den- 
sity variations as being due to (recently) exploded SNe and 
mar ginally trace their lig ht curve evolution. We note that 
also iMarchili et all l|2010h were able to identify and char- 
acterise SN 2008iz, within the nuclear regions of the nearby 
starburst galaxy M82, by carefully removing the diffuse nu- 
clear emission and fitting the resultant SN light curve. 

Although our flux density variation method is not very 
sensitive in the case of bright nucleus such as Arp299-A, 
where the diffuse emission dominates the overall radio flux 
density, it is well suited to be applied to sources like Arp 299- 
Bl, which at a distance of 44.8 Mpc has a luminosity around 
2.4 x 10 28 ergs -1 Hz -1 , thus being bright enough to provide 
a high RSN rate, but faint enough to enable their detection 
through flux density variations. 

The detection of new SNe in Arp299-A through this 
method (based on VLA measurements) is a challenging 
task: i) given the expected IR luminosity based SN rate, 
a fine sampling of the light curve is needed; ii) due to the 
brightness of the nucleus, new RSNe should have I/ pC ak > 
3.1 x 10 28 ergs -1 Hz -1 to become 5cr detections. 

Moreover, in the case of the A-nucleus, the expected 
typical RSN flux densities are comparable with the uncer- 
tainties in the nuclear flux density, making the SN detection 
impossible, unless very bright SNe (such as powerful type Iln 
or 1998bw-like events) occur. To overcome this situation, we 
have also tested the possibility of using high-resolution NIR 
AO observations for the detection of SNe within the inner- 
most nuclear regions of Arp 299-A. However, we were unable 
to detect any new SNe. Our estimated limiting magnitude 
for the SN detection at ~ 300 mas (or 70 pc projected dis- 
tance) from the centre of the nucleus A would allow the de- 
tection of typical CCSNe suffering from extinctions of Ay up 
to ~ 15 mag. At a larger offset from the nucleus (of 1 arcsec, 
or 200 pc projected distance) our observations were sensitive 
to typical CCSNe with Ay up to ~ 30 mag. 

The VLA epochs included in our sample were observed 
for different purposes and thus have different quality (i.e., 



sensitivity and resolution). Hence, a variability test that 
yields proper SN identifications through a correct light curve 
sampling has not been possible. The ideal case would be 
to have a specific observing program which provides homo- 
geneity in the observed frequency, resolution, and having a 
regular time span between observations of the same source. 
The EVLA represents a very good instrument to achieve 
these goals. However, its changing configuration represents 
a changing resolution through the year, and thus a finer sam- 
pling that allows the identification of type Ib/c SNe is again 
compromised. The upcoming eMERLIN, with a very high 
sensitivity and resolution, provides the best option to carry 
out this kind of studies. 

Our work has emphasized the importance of extinction- 
free measurements of the SN activity within the innermost 
nuclear regions of LIRGs. The SNe in the innermost nuclear 
regions of Arp 299 have been practically missed by the opti- 
cal and NIR studies due to the high extinction therein and 
the lack of sufficient angular resolution. We have shown that 
at radio frequencies, with relatively high-angular resolution, 
it is possible to detect SNe buried deep in the LIRG nu- 
cleus. Determining the complete CCSN rates in local LIRGs 
is of great importance also for interpreting the result of the 
CCSN searches at higher redshifts used to trace the SFRs 
as a function of redshift. 
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